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Abstract 

Motivated by multi-throat considerations, we study the phenomenological imphcations 
of a bulk axion in a slice of AdSs with a large extra dimension: k ~ 10~^eV, kR > 
0{1). In particular, we compare axion physics with a warped geometry to axions in 
flat compactiflcations. As in flat compactiflcation scenarios, we flnd that the mass of 
the axion can become independent from the underlying Peccei-Quinn scale. Surprisingly, 
we flnd that in warped extra dimensions the axion's invisibility, cosmological viability, 
and basic phenomenology remain essentially unaltered in comparison to axions in flat 
compactiflcations. 



1 Introduction 



The origin of CP violation within the standard model remains an essential theoretical 
problem. Current empirical evidence strongly favors the weak interactions as the source of 
observed CP violation through complex phases within the CKM matrix [1] . The lack of a 
Nambu-Goldstone boson interpretation for the r] meson, associated with the spontaneous 
breaking of a global U(l) flavour symmetry [2], suggests the existence of a non-trivial 
QCD vacuum structure [3] , which, in principle, violates CP. Extending the effective QCD 
Lagrangian to include the most general vacuum structure yields 

^eff = ^QCD + Q-^^F^^F^^a- (l) 

The effective strong CP violation parameter, 6, contains contributions from both QCD 
and quantum flavour dynamics, namely 

e = e + arg (det M) (2) 

where M denotes the quark mass matrix and labels the CP violation term arising 
solely from QCD. From the low energy effective field theory point of view, appears as 
a standard model input. Empirical measurements, such as those arising from limits on 
the neutron's electric dipole moment indicate [4], 

e < 10-^ (3) 

A priori, 6 could have a value anywhere on the interval [0, 27r] and the origin of this 
parameter's seemingly unnatural tiny value constitutes the strong CP problem (see [5-7] 
for reviews). 

Perhaps the most elegant solution to the strong CP problem rests on the Peccei-Quinn 
(PQ) mechanism [8]. In this scenario, becomes effectively promoted to a field - the 
axion - identified as a pseudo-Nambu-Goldstone mode associated with the spontaneous 
breaking of a global U(l)pQ symmetry [9]. As the axion dynamically relaxes to the 
origin of its potential, vanishes. Augmenting the QCD Lagrangian of eq.(l) by the PQ 
mechanism one finds, 

Ceff = >Cqcd + Id.ad^a + -^^^,F'""'F,ua (4) 
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where fpQ denotes the axion decay constant set by the scale of U(l)pQ breaking, and the 
parameter ^ arises as a model dependent factor. The axion obtains its mass [10] through 

QCD instanton effects that explicitly violate the U(l)pQ symmetry and thereby provide 
the axion with a potential, yielding 

ma ^ —r (5) 

JPQ 

where Aqcd ~ 250 MeV. Constraints from laboratory searches, astrophysics, and early 
universe cosmology bound fpq (see [7, 11] for reviews), 

10'GeV</pQ<10i2^eV (6) 



which, by eq.(5), implies. 



10-^eV<m„< lO-^eV. (7) 



We should emphasize that the PQ mechanism intimately hnks the scale of U(l)pQ break- 
ing, /pQ, and the axion mass. The PQ mechanism and the resulting axion present a 
solution to the strong CP problem, however, the particulars of the UV embedding re- 
main model dependent.^ 

One of the most striking axion scenarios results upon embedding the PQ framework 
within large extra dimensions [14, 15]. If the axion appears as a bulk pseudoscalar, novel 
features emerge [15], such as Kaluza-Klein (KK) mode interference, non-trivial axion mass 
relations, and increased energy dissipation among cosmological relic axions. With the ad- 
vent of warped compactifications [16], whereby large scale hierarchies can be understood 
through exponential red-shifting, several attempts were made in constructing warped 
axion models using both bottom-up approaches [17] and top down string constructions 
[18, 19]. In traditional warped extra-dimensional models that address the gauge hierarchy 
problem of the standard models, the first KK excitations of bulk fields occur near the 
electroweak scale, namely on the order of 10 TeV. However, generic flux compactifications 
allow for backgrounds with several "throats" [20] originating from a compact Calabi-Yau 
manifold. As demonstrated in [21], additional throats can provide new model building 
avenues for axions. In principle, the multiple throats could have widely different AdSs 



^ Cf. [12, 13] for the most well know implementations and see [6, 7] for reviews. 
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curvature radii and thus bulk fields would have their KK spectra dictated by the geom- 
etry of their throat. Motivated by these observations, we model-independently consider 
the extra-dimensional phenomenology of a bulk axion in a throat separate from the stan- 
dard model with an AdS inverse curvature radius of A; ~ 10~^ eV and with kR ~ 
For simplicity, we will assume that the standard model remains confined to a UV brane. 
In this scenario, the axionic KK excitations appear hierarchically lower than the elec- 
troweak scale. The setup provides a warped extension to the fiat extra-dimensional axion 
phenomenology studied in [15]. 



2 Axions in a warped background 

As a toy example and proof of principle, we will assume that the axion arises as a 
Nambu-Goldstone mode associated with a higher dimensional complex scalar field which 
spontaneously breaks a U(l)pQ symmetry. 



(8) 



" V2 

We will further assume that the vev remains constant over the entire AdSs slice. This 
scenario can arise through a simple multi-throat model where the standard model exists 
in its own throat (SM throat) while the complex scalar exists in a separate throat (PQ 
throat) as illustrated in figure 1. In the PQ throat we have, 

S^^M^J d'xy/\i\g^''dM^*dN^ (9) 

where upon decomposing ^ — r]{y) exp(ia(y)) we apply the boundary conditions 

daliR^O d,a\^y^O (10) 

which break the global U{1)pq symmetry and yield a constant solution (77) = v for the 
vev. The resulting axion, a{y), couples to the standard model via UV brane-localized 
interactions.^ 



^In [21], the UV brane localized interactions are mediated by exotic coloured fermions, Q,Q, with 



action Sq=J d'^xy^\gind\uv{^QLQR + ^*QrQl)- Under a U{1) chiral rotation of the Q fields, the 
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While generically we also expect a KK tower associated with the radial mode, we do 
not allow any radial mode- standard model interactions. Furthermore, the radial mode 
interactions appear through derivative couplings suppressed by powers of fpQ. As we will 
show later in this section, these couplings are further suppressed by a volume factor. As 
a result, for the remainder of this paper we will ignore the effects of the radial mode. 




(a) (b) 

Figure 1: Illustration of a two throat axion toy model, a) The standard model lives in left 
AdSs throat while a complex scalar field $ with phase a (the "axion") lives in the right 
throat. The two throats communicate via UV brane localized coloured fermions, Q,Q. 
b) Once the PQ symmetry is broken, the resulting bulk pseudoscalar axion effectively 
sees the standard model as UV confined matter. Each throat has its own inverse AdSs 
curvature denoted by ki and 

We stress that the above toy model only represents a proof of principle. We will take a 

model independent approach in that we will not consider any particular implementation 

of the PQ symmetry breaking mechanism within the AdSs set-up other than insisting on 

the constancy of the vev across the AdSs slice. In this paper, we only wish to explore the 

phenomenology of a bulk axion in a warped geometry. We will assume that the standard 

model remains confined to the UV brane and that the axion propagates in the bulk (as in 

figure 1 (b)). Since we assume that the axion potential only arises from QCD instanton 

axion coupling becomes transferee! to a gluon topological term of QCD on the UV brane, namely 
•Sint^ / d'^x {32tt^ fpQ)~^aFfj,i,F'^'^ which corresponds to the generic UV localized axion coupling we use 
in this article. 
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effects, we will omit bulk or boundary mass terms for the axion. Defining the line element 
within the AdSs space as, 



ds^ = e'''^^yr]t,^,dx''dx'' + dy^ = gMwdx^dx 



N 



(11) 



iVtii^ — diag(— 1, 1, 1, 1)) our action reads, 

1 



Seff — d xdy e 



-4ky 



9^^ a 



2 327r^ fpQ 



(12) 



The axion field, a{x,y), carries mass dimension [3/2] implying that [fpq] = 3/2. The five 
dimensional Newton's constant has been absorbed into the field definitions. Application 
of Neumann boundary conditions in the absence of bulk or boundary axion mass terms 
leads to the KK decomposition. 



a[x 



(13) 
(14) 



where 0o is a constant. Decomposing the action leads to, 
S= j d^xdy e-^^y 



(15) 



The kinetic terms. 



(16) 



contain no cross terms since the wave function profiles satisfy the orthogonality condition 

(17) 



Aky 2kyi (k = J) 



for all n,m - including the zero mode. We may write the wave function profiles as, 

'*[j,(!^e'=»)+rn(^e'=»)" 



0n 



Nr,. 



for n > with the normalization factor, 

tR 



(19) 
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The coefficients 63"^ ^-^e determined by, 



iSn) 
^2 



h\'{mn) 

and we now can re-express the 4-D action as 

n=0 n=l 



2J2{mn/k) + {m„/k)J2{mn/k) 
2Y2(mJk) + {mJk)Yi{mJk) 



(20) 
(21) 



(22) 



At this point, we will now make the definition, 

fpq = V27rRfp 



(23) 



leading to the parameter fpg which has canonical mass dimension 1. The parameter fpg 
implicitly depends on the 5-dimensional Newton's constant. Thus, a volume factor relates 
the two parameters in eq.(23). This allows us to set fpg to a hierarchically high scale. 

Like the flat extra dimensional counterpart [15], this construction solves the strong 
CP problem. Using the one-instanton dilute gas approximation, we have 



{F'^'^F^,) = -A^QCDsin ( -i 5^a,0„(O) + 9 ) 

\Jpq n=0 / 



which gives rise to the axion potential. 



V 



1 



+ 



9^ 



A 



n=l 



327,2 ^^QCD 



l-cos(-^5^a„0„(O) + e 

JPQ n=0 



The CP conserving minimum appears from. 



dV 

da„ 



2 , J. ^ 



^QCD sin ( -|- XI «n0n(O) + © 



(24) 



(25) 



(26) 



leading to the condition: 



/ ^ _ fpg - 9 

^"""^ ~ e 0(°Ho) 

(a„) = n > 0. 



pq n=0 



I e 2Z 



(27) 



We see that only the zero mode serves as the true axion. The higher dimensional axion 

mass matrix obtained from eq.(26) reads. 



da^dan 



2c2 A4 



QCD 



0m(O)0n(O) COS ( X anMO) + © 



327r2 f2 

Jpq 



pq n=0 



(28) 
(29) 



(a) 



and therefore we obtain, 



= rnlg ((/>„.(O)0„(O) + yl5mn) (30) 

with 

mpn = „^ — (31) 

yo = (32) 
yn = ^ n>0. (33) 

We can re-express 0*^'"^ (0) using appropriate Bessel function identities, 

An) _ Ji{mn/k) 



yielding 



^ \ k ) YAmJk) ^ \ k ) 



k / Yi{mn/k) V /c 
2.1 Axion KK mode mixing 



(35) 



The important difference between the usual 4-dimcnsional axion and its higher dimen- 
sional implementation centers on the effects of KK mode mixing given in eq.(30). It has 
been shown [15] that flat compactifications can lead to a divorcing of the axion zero mode 
mass from the PQ scale, leading to the approximate relation 



mo ^ min Qi? \ mpQ^ 



(36) 



for the lightest axionic state. The mixing of the axion KK modes, via the mass matrix 
eq.(30), leads to important differences between the physics of the usual four-dimensional 
axion and its higher-dimensional counterpart. In the standard four-dimensional axion 
scenario, the PQ scale not only sets the mass of the axion but also controls the strength 
of the axion coupling to the standard model. Observational constraints imply a PQ 
scale between 10^ GeV and 10^^ GeV corresponding to an axion mass window of 10~^ 
eV< mpQ < 10~^ eV. In the extra dimensional situation, where the axion appears as 
a bulk field, these standard arguments no longer apply: the would-be axion ao mixes 
with the KK tower via the mass matrix eq.(30). This mixing implies that mpQ no longer 



sets the lightest axion state. Using the property that the smallest eigenvalue of a real 
symmetric matrix is bounded above by the eigenvalues of any diagonal block, one can 
show (from the upper-left 2x2 diagonal block of eq.(30)) that the mass of the lightest 
axionic state is bounded above by both mpg and mi. Hence, rua < min(mpQ, mi), and we 
see that for nipQ > mi, the lightest axion mass decouples from mpg (and therefore /pq) 
- the lightest KK mass can set the scale. This result parallels the flat extra-dimensional 
case (see eq.(47) of [15]) where the inverse radius of the extra dimension can also set the 
lightest axion mass. 

Once the compactification radius and fpQ are chosen within the flat compactification 
scenario, the form of eq.(30) becomes fixed and leads to a definite KK mixing pattern. 
Warped compactifications lead to more freedom. In the large kR hmit {kR » 1), we 
have 



\^7r^ RrUn 

and, using asymptotic expansions for all the Bessel functions involved, we obtain 

2" 



(/)„(0) ^ -ne-^'''/^^/R^, 
The mass matrix eq.(30) becomes 



1 /rrinV' 
8 \ k ) 



(37) 



(38) 



/ 2irfefl 

-\/2,rafcil2m3e-'''=«/= 



^(mifl) + ■ 



'PQ 



•-PQ 



(39) 



This pattern substantially differs from the flat compactification scenario. The main new 
feature concerns the mass of the zero mode. Like the flat 5-dimensional case, the axion 
does not receive a mass term until mpq turns on. By contrast, in the highly warped case, 
the zero mode "axion" contributes a mass eigenvalue of 



mo ~ mpQVOjrkR. 



(40) 



In particular, this mass eigenvalue can be tuned relative to the flrst several KK states. 
Furthermore, the amount of mixing between the zero mode and the other KK states 
becomes tunable through the warping, kR. We note that if we choose kR <C 1 at flxed 
R, we recover the flat higher dimensional results. 
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We should also emphasize the nature of the mixing matrix at large KK number. In 
this case, 



and we find, 



Ar2 ^ esc^ ( ^ - ^ ) e'^^^ (42) 



which leads to 

(f)n^V2 (yfMe-''^^'^^ n > 1. (43) 
In this limit, eq.(30) now reads. 



= mpQ(7riJfc)e~''''" 



2 2 2 



'"PQ 

2 ^(n + l)^(Yfi)e-'^'°« 2 2 

2 2 2+ ^("+2)''('-^/«)e-'''" 2 

2 2 2 2 + ^("+3)''(Mfl)^~"'°-" 

""PQ 

(44) 



Apart from the exponential suppression factors, the mixing matrix approaches the fiat 
case form. We stress that the 0„ coefficients approach a constant as n — > oo. This feature 
will play an important role in the following sections. 



3 Decoherence, Invisibility, Energy Loss, and Detec- 
tion 

Perhaps the greatest potential threat to extra dimensional axion scenarios concerns the 
loss of invisibility. While each individual KK mode couples with 1/ fpQ (as in [15]), the 
"axion" is produced in a linear superposition of KK modes. The effective coupling to 
this superposition determines the axion's invisibility. Therefore, the effective interaction 
Lagrangian reads, 

jCi^t = l-LVN^Ff""'F.,a (45) 
where A'^ is a normalization factor defined by, 

N^J2(<^^nr (46) 

n=0 
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and 



a — 




(47) 



defines the "axion". We stress that only the zero mode state, a°, serves as the true 
axion in this set up as only of inherits the shift symmetry from the five dimensional 
theory. Since the entire linear superposition of the KK tower given by a' enters the 
efi^ective Lagrangian, the KK sum must be truncated - presumably at the cutoff of the 
5-D effective field theory. However, processes which produce on-shell axions should not 
include KK states above the characteristic production scale. 

If we consider axion processes at characteristic energy E, the effective field the- 
ory which results from integrating out axion mass eigenstates that exceed E reads, 



jCeff = K=oi9an? + T.2',;^^^M,nnaman + ^^K=t(t>nan. Ignoring the the effects 



of mixing from eq.(30), we can approximate rimax through run^^^ — E. We see that the 
gauge fields couple to the normahzed hnear superposition a' — {N{E))~^/'^T,'^^'^(j)nan 
with coupling strength ^jN{E)fpQ. The normalization factor N{E) = E^^q''^^ enhances 
the effective coupling. For clarity we have indicated the implicit energy dependence of 
N arising from the condition rrin^^ — E. Thus, the effective coupling strength of the 
linear superposition a' depends on the characteristic energy scale as does the composi- 
tion of a', and the determination of axion invisibility becomes much more subtle in an 
extra-dimensional scenario. Furthermore, the state a' produced in interactions with vis- 
ible matter is manifestly not a mass eigenstate, and so neither its mass nor its lifetime 
are well-defined. 

Since N depends on KK number, and since the kinematics of a given process fix the 
number of modes included in a', the effective coupling in eq.(45) grows with energy. In 
general, eq.(30) creates a misalignment between mass and interaction eigenstates such 
that 



where a denotes the mass eigenstate and U diagonalizcs the axion mass matrix, cq.(30). 
From the 4-dimensional perspective, at a given energy scale, a' appears as a field in the 
interaction basis. As the state a' propagates, the individual KK modes interfere, creating 
an axion oscillation phenomenon. Again, we stress that a' consists of a superposition of 




(48) 



n=0 
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mass eigenstates and does not have a well-defined lifetime or mass. We will examine two 
important physical consequences: energy loss and direct detection. 

The strongest bounds on the 4-dimensional Peccei-Quinn scale arise from stellar cool- 
ing constraints. Weakly coupled axions provide new channels for stellar energy loss and 
the total luminosity in axions must not upset stellar evolution. The current lower bound 
on the Peccei-Quinn scale in the usual 4-dimensional axion scenario from HB stellar 
cooling reads (see [11] and references therein), 

/pQ>2xl09GeV. (49) 

This bound applies provided that the axion mass does not wildly exceed the internal 
temperature of the star and the coupling remains sufficiently weak [11]. In the higher di- 
mensional analogue, the effective coupling contains energy dependence from the inclusion 
of KK modes. Thus, eq.(49) becomes modified, leading to 

j^^ > 2 X 10^ GeV (50) 
V N{E) 

where N{E) denotes the normalization factor of eq.(46) which includes modes up to Umax- 
For simplicity, we will assume that the mode superposition within a' contains states up 
to a cutoff defined through the internal stellar temperature via mn^^^ — T. We will also 
make the further approximation that the usual 4-dimensional kinematical phase space 
constraints apply for a' production.^ 

In the warped case, we can provide an estimate on the energy dependence of N{E). 
If we approximate the mass spectrum by the KK masses (again, ignoring the effect of 
eq.(30) on the mass eigenvalues), at large KK number we have 

E^m^^^^nke--^'' (51) 

which leads to 



"^Strictly speaking, wo expect that thermal corrections to the phase space will play a role in the 
production of the more massive states in the superposition, as shown in [22]. Since we crudely truncate 
the tower at T, for the purposes of this paper, we ignore these effects. 
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Prom eq.(46), and using the asymptotic approximation for all eq.(43), we can obtain 
an approximate explicit formulation oi N{E) in terms of the kinematic cutoff, 




N{E) 



Thus, at a fixed energy, the effective Peccei-Quinn scale for a' becomes 

fpQ 



(53) 



eff 

PQ 



(54) 



E V2RE 

where we assume RE 1. In order to satisfy HB stellar evolution constraints (assuming 
core temperatures of ~ 8keV [11]), we require 



JPQ 

implying. 



£=8keV ^ 



> 2 X 10^ GeV (55) 



> 2 X 10^ GeV. (56) 



V(2i?) -SkeV 

Note the above result obtained for a warped compactification gives comparable bounds 
for flat compactiflcation [15,22]: 

rimax = EE (57) 

and thus, 

AT^ = 1 + ^ 2 fti 2ER. (58) 

n=l 

Both the flat and warped compactiflcations give the same normalization factor at flxed 
energy. 

While the enhanced coupling must satisfy the astrophysical energy loss constraints, the 
effects of the superposition within the a' fleld lead to phenomenological consequences in 
direct detection searches. Again, the relevant quantity is not the individual KK modes nor 
the suppressed coupling fpg, but the superposition a' with coupling /p'g . Even though 
the extra-dimensional scenario leads to an enhanced coupling, the KK modes interfere 
during propagation. Only the linear superposition a' couples to standard model flelds. 
It has been shown [15] in flat compactiflcations that destructive interference is crucial in 
re-establishing invisibility as it reduces the expected flux measured from distant sources. 
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In order to determine what fraction of the flux produced by a distant source can be 
detected, we must calculate the survival probabihty of the a' field itself. The amplitude 
for individual KK axion transitions is given by, 

Ak^i{t) = J2UiiU:kexp (59) 

where denotes the mass eigenvalue arising from eq.(30), and U labels the unitary 
matrix that diagonalizes the axion mass matrix. We have ignored the possibility of axion 
decay throughout. The probability that a' remains preserved during propagation reads, 

2 



y^^(t)l(t)kAk^l{t) 



(60) 



k,l 

Therefore, the expected measurable flux arriving from a distant source appears as. 



$ = (P„'_„')*o (61) 

where $o denotes the flux as calculated using the enhanced coupling f^Q in the absence 
of decoherence. Since {Pa'-^a') < 1 we see that interference serves to reduce the the total 
measurable flux. 

As an estimate of the decoherence time we consider the time it takes for the largest 
mass eigenstate in the linear superposition, a', to step completely out-of-phase with the 
zero mode. Applying this out-of-phase condition to eq.(59) yields, 

= cos f!^i-l^tol (62) 



2p 



and we flnd that. 



* __£Q_ 2nkR (no\ 
~ {uT^kf' ■ ^^^^ 

The parameter Tq, defined through Tq = to{m%Q/2p) (where p denotes the usual 3- 
momentum), gives a dimensionless time element. In the fiat case limit, eq.(62) yields, 

2 

(64) 



in agreement with [15]. 

Note that Tq denotes a "decoherence time" defined by Pa'^a'(To) ~ {Pa'^a') ^ i-e. 
the time taken for Pa'^a'{t) to reach its time independent average value. The scaling 
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argument proceeds subtly. We know that the axion remains in the a' state as mpQ — > 
since the mass matrix of eq.(30) becomes diagonal in this limit and, as a result, no mode 
mixing can occur (i.e. the transition amplitudes become diagonal, Ak-^i = 6ki). This 
leads to the complete lack of decoherence and reproduces the usual 4-dimensional axion 
scenario. As mpg —> we find that Tq — > 0, implying a vanishing decoherence time. This 
result remains consistent since Pa'-^a'{T) ~ {Pa'-^a') — > 1 in this limit. Thus, even though 
Tq ^ as mpQ ^ we also have {Pa'^a') 1 and hence decoherence does not occur. ^ 
The largest KK mass in both the fiat and warped cases govern eq.(63) and eq.(64). 
We may re-write the decoherence time as 



which holds for both cases. Approximating the energy of the linear superposition through 
the mass of the heaviest mode we find, 

which is independent of the details of the compactification. Thus, provided that many 
modes appear in the superposition, large warped and flat compactifications lead to the 
same decoherence time at fixed energy. 

Furthermore, we can obtain a semi-analytic expression for the resulting time indepen- 
dent average value (Pa/ _>„'). In the flat compactification at fixed mpg and fixed R, it has 
been shown [15] that {Pa'-^a') ~ ^/n. This behaviour follows from expanding eq.(59) and 
eq.(60). Applying these results to the flat compactification scenario gives 



n 



a ^ -1-26 cos 

i i<j 



[mf - mf\t 



(67) 



2p 

The 0{1) coefficients a and h result from from the unitary matrix that diagonalizes the 
flat extra-dimensional axion mass matrix. At i = eq.(67) gives the expected result 



Pa'-.a' ^[\]n'^l. (68) 



*Wc stress that our definition of decolierence differs from [15] where the authors define deeolierenee 
by fixing Pa'^a'iTa) = 0.90 for all values of mpg and R. Their definition requires a non-trivial mode 
number renormalization in which the energy dependence of fpQ , of the a' composition - and therefore 
the energy dependence of Pa'^a' ~ is not manifest. 
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At later times, the cosine terms no longer add coherently and thus eq.(67) reduces to 



The warped case proceeds in a similar fashion. By observing that eq.(44) has the 
same structure as the resulting flat compactiflcation matrix (apart from an overall mul- 
tiplicative factor), we see that eq.(67) applies. In particular, at large mode number, the 
essential difference between the warped and flat compactiflcations centers on the density 
of KK states. At a fixed energy, the warped compactiflcation contains ~ TTkRexp{7ikR) 
more modes than the flat case. Eq.(67) then tells us the scaling in the warped case 
appears as. 



indicating an exponentially larger value for the expectation of Pa'^a' at fixed n as com- 
pared to the flat compactiflcation result. However, if we compare the results eq.(69) and 
eq.(70) not at flxed n, but at flxed energy we flnd that in both cases 



since the warped case yields the approximate relation E ~ nnkR exp{T: kR). Remarkably, 
at flxed energy, both the warped and flat extra dimensional axion scenarios lead to a 
similar value for {Pa'-,a')- Thus, while it appears that at flxed mode number warped 
compactiflcations permit decoherence times that are substantially longer than those oc- 
curring in the flat case limit, physical processes cutoff the mode number in both cases to 
yield a result that is largely independent of the compactiflcation. 

As examples, we show the effect of warping on the decoherence time relative to flat 
compactiflcations in flgures 2 and 3. In both cases, we have taken R — 1/(10~^ eV) 10 //m. 



In figure 2, where a fixed number of 30 KK modes are used, it appears as though 
the reduced mixing in the warped compactification scenario (see eq.(30)) reduces the 
axion's ability to decohere. Figure 2(a) displays the effect with a mild warp factor of 
kR — 0.5. In this case, we find that the flat compactiflcation large time limit yields 
{Pa'^a') ~ 3.3 X 10~^ while the warped compactiflcation gives a similar value of {Pa'^a') ~ 
3.9 X 10~^. In figure 2(b) we use a larger warp factor of kR = 1. The fiat compactification 




(69) 




(70) 



fpQ = W GeV. 
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reproduces {Pa'->a') ~ 3.3 x 10~^ in the large time limit. On the other hand, the warped 
compactification large time limit yields {Pa'^a') ~ 0.23, marking a significant departure 
from the flat compactification result at fixed mode number. 




Figure 2: The axion survival probability with kR = 0.5 and kR = 1 with 30 KK modes 
included. Plots a) and b) displays the flat compactification (dashed-dotted) along with 
the warped compactification (solid line). In all plots R = lO^eV"^, and fpq = 10^^ GeV. 

However, in physical processes, the maximum mode number within the linear super- 
position is fixed by kinematics. While the axion mass matrix in the warped case appears 
with suppressed mixing relative to the fiat case, the KK spectrum also appears more 
dense. The increase in mode number in the warped case makes up for the suppressed 
mixing leading to eq.(71) - i.e. the same result as derived in the fiat compactification 
scenario. Figure 3 displays the results at fixed energy, again for kR = 0.5 and kR = 1. In 
this case we have taken the mode number cutoff determined through E = 0.1 eV in both 
cases. We see that Tq and the final value for {Pa'~*a') appear approximately the same. 

4 Comment on Cosmological Relic Axions 

Cosmological relic axions provide strong constraints on the invisibility of the axion in 
the usual 4-dimensional PQ mechanism. During the universe's early thermal history, 
the universe passed through the QCD phase transition, at time tqcD, at which point 
QCD instanton effects established an axion potential where none existed previously. By 
naturalness, we expect that the axion would find itself displaced from the minimum of 
its potential by an (9(1) fraction of fpq. At this point, the axion would begin to oscillate 
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Figure 3: The axion survival probability with kR = 0.5 and kR = 1 at a fixed energy 
of 0.1 eV. Plots a) and b) displays the fiat compactification (dashed-dotted) along with 
the warped compactification (solid line). In all plots R = lO^eV"^, and fpq = 10^^ GeV. 
Note that Tq ^ [rripq/ E)"^ and {Pa'^a') ~ l/ER in both cases. 

about the minimum via 

d?a .da r, , 

— + 3H{t)— + mla = 0. t > Iqcd (72) 

with the initial condition aitQco) ~ fpQ- These relic oscillations, while damped through 
Hubble expansion, continue to store energy and must not exceed the present day critical 
density. Cosmological considerations provide a constraint from above on fpq in the usual 
4-dimensional PQ mechanism, namely fpq < 10^^ GeV. 

As displayed in figure 4, the extra dimensional situation proceeds differently. As each 
KK mode has a potential prior to the QCD phase transition and following [15], we can 
approximate each KK mode as sitting at the minimum of its potential at tqco- By 
contrast, the zero mode (which represents the true axion) does not receive a potential 
until the instanton effects turn on at tgcD- Thus, we expect at tgcD that the zero mode 
will have a natural displacement from its minimum, i.e. ao(tQCD) ~ /pq, o^iitqcD) = 
for / > 0, and dai/dt = for all /. In the extra dimensional case, eq.(72) reads, 

d'^di , ^ do,] , ,9 /_ X 

^ + 3i/(t)^ + M2afc = t>tqcD (73) 

where denotes the non-diagonal axion mass matrix. Once the zero mode begins 
to oscillate, the non-diagonal mass matrix in eq.(73) induces sympathetic KK mode 
oscillations that potentially have the ability to alter the cosmological constraints obtained 
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Figure 4: Higher dimensional induced KK oscillation plot as a function of time (dimen- 
sionless) with tqcD = 1; fn-pq = 4 x 10"^ eV, R = l/{2 x 10~^eV). Plot a) shows the 
first three KK mode oscillation with flat compactification (dashed-dotted) lines. The 
solid line denotes the usual 4-dimensional axion scenario. Plot b) displays the effect of 
warping with kR = 2 among the first three modes. Again, the solid line denotes the 
usual 4-dimensional axion scenario. 

in the usual 4-dimensional case. It has been shown in [15] that for certain parameter 
choices the oscillating KK modes can dissipate the oscillation energy more quickly relative 
to the 4-dimensional case, allowing fpq to become as large as ~ 10^^ GeV. However, once 
constraints from low energy phenomenology and gravity are applied, it can be shown [15] 
that the extra-dimensional situation dissipates the oscillation energy at the same rate as 
the usual 4-dimensional scenario. Remarkably, cosmological constraints on the flat the 
extra-dimensional axion scenario lead to the same viability as the usual 4-dimensional 
axion. 

Following [15], the results of the fiat extra-dimensional case can be readily extended 
to warped compactifications. We can transform to the mass eigenbasis such that eq.(73) 
appears as a set of uncoupled ordinary differential equations, 

d'^ai 3 da, 2~ n f7A\ 
— r^r H ; — h mfa? = 0, 74 

where we have assumed a radiation dominated epoch, H{t) = 3/2t, and we have de- 
fined di = di/ fpQ. For each uncoupled differential equation in eq.(74), we can use the 
dimensionless time element r = rriit leading to, 

d^di 3 ddi 

^ + 27^7 + "' = °' 
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The initial conditions now read, 

rin 7 

= (76) 



~ I \ TT ^^l- 



T=TO 



where U diagonahzes the axion mass matrix, eq.(30). Using the initial conditions in 
eq.(76), the general solution to eq.(75) reads, 

ai{T) ^ -^Uoi4^'T-'/'j{ro;r) (77) 

where 

j{ro; t) = J-5/4(to) Ji/4(t) + J5/4(to) J-i/4(t). (78) 

This solution matches the results given in [15] up to the unitary matrix element, Uqi. 
Defining the dimensionless mass eigenvalue ihi = mi/mpq the total energy density can 
be written as, 

with p = pI {m^pQfpo)- More succinctly, by defining t — mpgi, the total energy density 
becomes. 

Pit) = -^ioi''^^ ^o>f {jimio; rhiif + 3\mk\ mif) . (80) 
In the large time limit ihii » 1 eq.(80) becomes, 

p{i) = lx{i,)V'jn-'/' (81) 
where the time independent coefiicient Xito) reads 

^(^~o) = U^irhf {[J5/4{rhiio)]^ + [J-5/4{rhiio)f + V2J5/4(mzto)^-5/4(m/to)) • (82) 

Again, this result appears the same as in the fiat compactification [15] up to the unitary 
matrix element. Since the energy density for the four dimensional case in the large time 
limit has the same form as eq.(81) except with 

X4D{to) = [J5/4(to)]' + [J-5/A{io)? + ^275/4(^0)^-6/4(^0) (83) 

the ratio of the energy densities can be expressed as [15], 

r ^ ^ = (84) 
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As in the flat case, can significantly deviate from unity leading to a weakening of 
the upper bound on fpg. This deviation occurs only in the limit where iorhi < C(l) 
(i.e. the Bessel functions in eq.(82) do not approximate their asymptotic limit) for a 
significant number of modes which also includes the mass eigenvalue contributed by the 
mode Qq. Nevertheless, saturates at unity which implies that extra-dimensional axions 
- warped or flat - remain viable cosmologically. In figure 5 we have displayed three cases, 
to = 0.1,1,10, with mpqR = 1. We find that as the warping becomes increased, the 
deviation of Vp from unity increases as well, provided that the mass eigenvalue contributed 
by oo satisfies torn < 0{1). While the plots in figure 5 demonstrate the potential ability 
of warping to further increase the viability of the extra dimensional axion scenario, we 
caution that the large departure of from unity crucially hinges on certain parameter 
choices. In the fiat case, the condition required for < 1 reads Iqcd/R ^ - 
a condition requiring R ~ 10^°eV~^! This result wildly contradicts constraints from 
gravity. The warped case proceeds more subtly. If the zero-mode in eq.(30) contributes 
a typical mass eigenvalue of ~ exp{—7rkR) to the KK spectrum then, provided that 
tgcokTr exp{—7rkR) < 0{1), a departure of from unity can be expected. However, in 
the large warping limit {kR > 1), the zero- mode clearly contributes ~ mpQV^nkR (see 
eq.(30)) to the KK spectrum thereby requiring tQcnif^PQ^'^T^kR) < 0{1) in order for 
us to expect Vp < 1. Using Iqcd ~ 10~^sec ^ 10^°eV~^, and kR ~ 1 with R ~ 10 /xm, 
we find no significant deviation of from unity unless mpg < 10~^° eV (which generally 
does not include an enhancement factor larger than 0{1)). Thus, over the parameter 
range of interest, we can expect that warped extra dimensional axions dissipate the 
relic oscillation energy no more quickly that the usual four dimensional scenario. This 
result matches the flat extra dimensional case. We stress the remarkable central point: 
the extra dimensional axion scenario - whether warped or flat - continues to remain 
cosmologically viable, dissipating the rehc oscillation energy at at rate no less than the 
usual four dimensional axion scenario. 
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Figure 5: Ratio of the energy densities for to = 0.1,1,10 as a function of the warp factor 
kR with rripQ = 1/R. The k = hmit yields the flat case results in [15]. 

5 Conclusions 

In this paper, we have compared the basic phenomenology of axions embedded in a flat 
extra dimension with axions embedded in a warped extra dimension. Given the recent 
interest in multi-throat scenarios [20], large warped extra dimensions may provide new 
model building avenues. Furthermore, current and planned experimental axion searches 
will probe new parameter ranges in axion physics [23]. 

In a flat extra dimension, axion oscillations cause decoherence which leads to invis- 
ibility. We have found that warping the extra dimension controls the axion oscillations 
such that the decoherence length becomes tunable in KK number. While this result 
appears to give greater freedom in establishing the level of axion invisibility, kinematics 
govern the number of modes that must be included in the axion linear superposition. 
As the warped case demands that we include more modes relative to the flat case at a 
given parameter choice, we surprisingly find that the overall effect leaves the invisibility 
unchanged relative to flat compactifications. 

Furthermore, in both warped and flat extra dimensions, kinematics also govern the 
effective coupling, fp^g. While each mode couples with safe ~ (pn/ fpQ, the on-shell state 
couples to the standard model with the much larger fpg ~ fpq/ a/ N{E), where N{E) 
is a function of the production energy. Thus, the production kinematics play a central 
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role in determining axion's coupling to the standard model - for both warped and flat 
compactifications. Constraints derived from energy loss mechanisms must take this effect 
into account - it is not the coupling of an individual mode the establishes the bound 
but rather the coupling of the superposition. Within extra-dimensional scenarios, the 
energy dependent effective coupling, fpg, controls astrophysical axion production, and 
therefore controls the source's axion luminosity. On the other hand, decoherence, arising 
from KK axion oscillations, sets the expected measurable flux. In the warped geometry, 
the AdS curvature, k, plays a crucial role in both effects and yet for processes that involve 
large mode numbers, the expected phenomenology remains essentially unaltered from 
flat compactifications. It would be interesting to perform a detailed analysis involving 
constraints from astrophysical sources on axions in the warped extra dimension scenario 
we consider. 

We have also found that warped compactifications remain cosmologically viable. As 
in the flat case, we have found that the sympathetic KK mode oscillations induced by 
the displaced zero-mode dissipate the axion oscillation energy at least the same rate 
as the usual four-dimensional axion over most of the parameter range. If we include a 
large warp factor {kR > 1) and take mpQiqcD ^ we find the possibility that the 

KK modes dissipate the relic oscillation energy more quickly than the four-dimensional 
axion scenario. While this result in principle allows a larger value of fpq relative to the 
usual four-dimensional bound (~ 10^^ GeV), the resulting factor by which the dissipation 
becomes increased is generally insufficient to allow a significant increase of fpq. 
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